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Abstract. We study convective motions taken from hydrodynamic simulations of rotating proto-neutron stars 
(PNSs) with respect to their ability to excite a dynamo instability which may be responsible for the giant neu- 
tron star magnetic fields. Since it is impossible to simulate the magnetic field evolution employing the actual 
magnetic Reynolds numbers (Rm) resulting from the hydrodynamic simulations, (smallest) critical Rms and the 
corresponding field geometries are derived on the kinematic level by rescaling the velocity amplitudes. It turns out 
that the actual values of Rm are by many orders of magnitude larger than the critical values found. A dynamo 
might therefore start to act vigorously very soon after the onset of convection. But as in general dynamo growth 
rates are non-monotonous functions of Rm the later fate of the magnetic field is uncertain. Hence, no reliable 
statements on the existence and efficiency of PNS dynamos can be drawn without considering the interplay of 
magnetic field and convection from the beginning. Likewise, in so far as convection inside the PNS is regarded 
to be essential in re-launching the supernova explosion, a revision of its role in this respect could turn out to be 
necessary. 
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1. Introduction 

The origin of the neutron star (NS) magnetic fields 
is subject of scientific debate since the discovery that 
\ the pulsar mechanism is based on the existence of the 
strongest magnetic fields seen in the universe. Though 
the idea of magnetic flux conservation during the col- 
lapse of the NS progenitor came up quite naturally, 
doubts as to whether that mechanism explains the ob- 
served fie ld strengths were raised alm ost simultaneously 
(see, e.g.. lThompson fe Duncanlll993|h That macroscopic 
currents - as almost everywhere else in the universe - 
are the cause of the NSs' magnetic fields seemed obvious; 
thus attention has been dra wn to the idea of a dynamo. 
iFlowers fe Rudermanl l|l977l) first pointed out that such an 
instability can start to act immediately after the neutron 
star's birth, and last until its internal motions disappear. 

The increa sing understanding of the supernova mech- 
anism (e.g., iBurrows fe Frvxelil Il992l and references 
therein) and the accompanying process of NS creation 
strengthened the idea of a convective proto-neutron star 
(PNS) phase. This very first epoch in a NS's life starts 
immediately after the collapse when negative radial gradi- 
ents of both entropy and lepton number are created, thus 
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in gen eral enabling Ledoux conv e ction (see, e.g.. lEpsteinl 
(1979) and Burrows fc Lattimerl Jl 986ii and lasts about 
half a minute. iThomnson fc Duncan 1)19931 ) sketched the 
possibility of dynamo action caused by turbulent convec- 
tive motions in the PNS. They claimed that as these vig- 
orous motions combine with differential rotation, a mean- 
field dynamo of a— ft type may act in regions where the 
rotation is capable of influencing the convection signifi- 
cantly. By use of a mixing-length approach they estimated 
the convective velocity to be about 10 8 cm s^ 1 and con- 
cluded from equipartition that such a dynamo may am- 
plify a seed field up to 10 15 G during the first 30 seconds of 
a NS's life. In contrast, small-scale dynamo action with- 
out generation of a coherent large-scale field should be 
expected where the influence of the rotation on the con- 
vection is negligible. 

In two complement ary papers, iBonanno et al.1 l|2003h 
an d llTrnin fc Oil l)2004h recently presented studies on field 
generation in PNSs where the former deals with mean- 
field dynamos and the latter with small-scale ones. 

In their strongly simplified model, IBonanno et alJ 
(2003) assume the mean e.m.f. to be constituted by 
isotropic a and (5 effects only and adopt purely 
radial dependences for differential rotation, the a— 
parameter and the turbulent magnetic diffusivity. Like 
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Thompson fe Duncan! l|l993j) and based on the results of 
Miralles et alJ l)2002|b thev argue that a is negligible in an 
inner, entropy-gradient driven convective region, but pro- 
portional to the rotation rate £1 in an overlaying, lepton- 
gradient driven convective shell. They derive critical val- 
ues of fl above which dynamo action in either the a 2 or 
the a—il regime is possible and conclude that the vast ma- 
jority of PNSs will show mean-field dynamo action. 

The work of lUrpfn fc Giil l|2004h is focused on small- 
scale dynamos which are supposed to emerge from both 
convective regions. As the magnetic Reynolds numbers are 
huge even when calculated with the scales of the convec- 
tive eddies, at least a necessary condition for this supposi- 
tion is surely satisfied. To get estimates for the magnetic 
field strength at the end of the convective phase, equipar- 
tition is employed thus arriving at magnetic field strengths 
of rs 10 13 G in both convective regions for the largest tur- 
bulent scales Ri 1 ... 3 km. 

Both studies state that given the estimated convection 
amplitudes, dynamos can act because the velocities are 
overcr itical with respect to these dynamos. lUrpm fc Gil 
(2004) do not estimate the real level of 'overcriticality', 
i.e. the 'distance' between the real a nd the critical velocity 
(e.g. as a ratio of amplitudes) at all. lBonanno et all l)2003|) 
provide a hint that in terms of the mean-field coefficient 
a the real flow might be overcritical by a factor up to 1000 
but refrain from discussing the consequences on the basis 
of the behavior of growth rates. 

But, as in general the growth rates do not simply grow 
with growing overcriticality (in the above sense) no reli- 
able conclusion on dynamo action can be drawn without 
calculating them based on realistic velocity amplitudes. 
Moreover, it appears desirable to overcome the drawback 
of using ad-hoc assumptions on the convection and in- 
stead to make use of 'realistic' convective velocity pat- 
terns obtained from hydrodynamic simulations of the PNS 
stage. By feeding them into dynamo calculations 'as they 
are', a better understanding of the different PNS layers' 
contributions to the field generation and hence the rela- 
tive importance of small and large-scale dynamos could 
be gained. 

Along with the problem of magnetic field generation, 
the comparatively short period of convection in PNSs has 
been taken to play an important role in the puzzle of an- 
other major question: How can a sufficiently large amount 
of heat and / or neutrinos be transported into the region be- 
hind the shock front to enable the re-energized supernova 
explosion? 

Recently, doubts have been reinforced about whether 
the PNS convection alone is sufficient to accomplish this 
and new hints have been provided that convection in a 
second outer sp herical shell (50 . . . 100 km) could be cru- 
cial instead fsee lJanka et al.l2004l and references therein). 
1 But at least in the case when the thermal energy stored 

It seems conceivable that a magnetic field is being gener- 
ated in this outer zone, too, parts of which could later consti- 
tute the crustal magnetic field by virtue of fallback accretion. 



outside the neutrinosphcre and transported by this con- 
vection towards the shock is not sufficient to re-launch 
it, the transport of neutrinos trapped inside the neutri- 
nosphere by virtue of the coupled action of the inner 
(PNS) and outer convective zones could help: If the over- 
shooting zone of the former and the undershooting zone of 
the latter overlap, a sufficient number of neutrinos would 
be enabled to escape from the PNS and to reach the post- 
shock heating zone. 

Up to now, all attempts to cope with the problem of 
the convection-driven explosion were made with the influ- 
ence of the magnetic field neglected. But, if success or fail- 
ure depends sensitively on details of the convective phase, 
a clear understanding of how early and how fast the mag- 
netic field starts to act back on the motions which gen- 
erated it (perhaps in bo t h con vection zones) seems to be 
essential. iMiralles et alJ l)2002|) were the first to include 
the influence of an imposed magnetic field on the condi- 
tions for the onset of the relevant convective instabilities 
and found that only extremely strong (10 15 G) magnetic 
fields suppress the occurrence of convection completely. 

Since the mid nineties fully nonlinear hydrodynamic 
supernova simulations also resolving the PNS started to 
provide m ore detailed insight into the convection there 
(see, e.g., lKriilll 997t I.Tanka fc Keill ll 99St iFrver fc Hegel 
Sffll I.Ta^kalEool I.Tanka et a,lJl20f)l. 

2. Kinematic dynamo action of PNS convection 

The purpose of this paper is to study the dynamo-related 
properties of velocity profiles found in hydrodynamic PNS 
simulations. To start with, we restrict ourselves to the 
kinematic approach, that is, we solve the induction equa- 
tion with prescribed velocities v. It reads in normalized 
form for a homogeneous medium 

B = A B + Rm curl(t) x B) , divB = (1) 

where Rm = UR/rj is the magnetic Reynolds number with 
a characteristic velocity U, the magnetic diffusivity 77 and 
the radius R of the PN S. We define v by sna pshots 2 taken 
from the simulations of lKeill l|l997l) (see also l.Tanka fc Keil 
1998) who performed 2D (axisymmetric) simulations of 
rotating PNS convection with a model including a com- 
pletely radial neutrino transport scheme and employing 
the EOS of Lattimer and Swesty without accretion. To 
check both an early stage possibly relevant for the explo- 
sion and a stage with fully developed convection perhaps 
determining the later NS field, we took the velocity pro- 
files at 30 ms and 0.9 s after bounce. 

The geometrical structures of the convective velocities 
v are shown in Fig. ^ Table ^ summarizes their main 
qualitative features. At the late stage the azimuthal ve- 
locity component (i.e., the differential rotation) is in a 
large portion of the PNS's volume nearly independent of 
z, as a consequence of the Taylor-Proudman theorem. The 

2 Using the originally time-changing convection instead of 
snapshots would not influence our conclusions qualitatively. 
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Fig. 1. Convective velocities simulated bv lKeill for 30 ms (left) and 0.9 s (right) after bounce in cm s 1 . In the 

left semicircles the isolines of the azimuthal velocity and in the right ones the streamlines of the meridional velocity 
are shown. Note that the former is measured with respect to a frame in which the flow has zero angular momentum. 
In the right semicircles background colors provide the modulus of the meridional velocity. 



Table 1. Qualitative features of the convective veloci- 
ties: R - PNS radius, R l7 R - estimated inner and outer 
radii of the convective shell; Ng, N r - number of prominent 
convective cells in meridional and radial directions; w max - 
maximum convective, t)™ ld - maximum meridional veloc- 
ity. Rm - magnetic Reynolds number based on the r.m.s. 
value of the velocity with respect to the full PNS volume 
and an electric conductivity of 10 24 s _1 . 
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early velocity pattern is dominated by only a few convec- 
tive cells per hemisphere whereas the late one shows more 
and, hence, smaller cells. 

As the velocity patterns are reficctionally symmetric 
with respect to the equatorial plane and symmetric about 
the PNS's axis of rotation, any eigensolution of is 
either equatorially symmetric (S) or antisymmetric (A) 
and its spherical components vary like cos my or smmtp, 
in = 0, 1,2, . . . where tp is the azimuth with respect to 
the rotation axis. For non-decaying solutions the case 
m = (axisymmetry) is excluded by Cowling's theorem. 
As a consequence, if dipolar dynamo solutions, i.e. in our 
case, Sl-solutions, occur their axes are bound to lie in 
the equatorial plane. Note that this restriction is only due 
to the axisymmetry of v and does not apply to a realis- 



tic situation with 3D convection. Then, it is well possible 
that a dominating axisymmetric (aligned) dipolar field, 
accompanied by non— axisymmetric constituents is gener- 
ated. We restrict ourselves here to the case m = 1 only 
as from experience we expect those modes to be the most 
easily excitable ones. Both equatorial symmetries are in- 
cluded, that is, for either velocity the Al and SI solutions 
are considered. 

For the boundary conditions at the surface we assumed 
the surroundings of the PNS to exhibit a significantly 
lower electric conductivity than its interior. We modeled 
this condition by requiring that the magnetic field is at 
the surface continued into a vacuum field. 

As we cannot hope to solve Eq. with the velocity 
snapshots' real values of Rm (see Tabled we restrict our- 
selves instead to determining the so-called critical values 
of Rm, Rm cr it, for which the total magnetic field energy 
(or a suitable temporal average of it) neither grows nor 
decays. The corresponding magnetic field solution is then 
called marginal. To accomplish this, the magnitude of the 
velocity was considered freely rescalable while its geom- 
etry remained fixed. That is, we simply multiplied the 
complete velocity field with a factor which was adjusted 
such that a marginal solution emerged. In mathematical 
terms, we sought solutions of Eq. Q with an exponen- 
tial ansatz for the time dependence of B, i.e., d/dt — A, 
where the real part of the time increment A has to van- 
ish. In general, this is not equivalent with stationarity of 
the marginal solutions as a non- vanishing imaginary part 
of A may occur causing oscillating marginal solutions. In 
our situation, however, in which we excluded axisymmet- 
ric (rn = 0) solutions, the oscillations introduced in this 
way are restricted to uniform rotations of the field pattern 
about the axis of rotation, where the imaginary part 3{A} 
defines the rotation rate. Hence, the magnetic field energy 
is for any marginal solution a constant quantity. 
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Table 2. Critical magnetic Reynolds numbers and char- 
acteristics of the marginal eigensolutions for the velocities 
considered. E po \ and E tm are the energies of the poloidal 
and toroidal parts of the field, respectively, where the 
poloidal part comprises the externally observable multi- 
pole constituents. The toroidal field is zero outside the 
PNS and its field lines lie completely in spherical surfaces. 
(l, n )max are the maximum spherical harmonic order and 
number of radial zeros, respectively, included in the nu- 
merical model; (Z,n)dom correspond to the scales domi- 
nating in the energy spectra of the solutions. 
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Our actual technique to determine i?m cr j t together 
with the corresponding marginal solution was to integrate 
Q as an initial value problem. We started with values of 
B specified according to the desired symmetry type (Al 
or SI) of the solution, but arbitrary otherwise. In every 
time step we adjusted Rm (automatically) such that the 
total magnetic energy approached a stationary value. As 
no marginal solutions exist for which the total magnetic 
energy is oscillating no modes of a given symmetry type 
can be overlooked the i?m cr it of which could possibly be 
smaller than the detected one. 

We used a spectral code based on the magnetic field's 
free decay modes, specified for vacuum boundary condi- 
tions. Each of these modes is characterized by a spherical 
harmonic F ; m ( , (9, <p) providing the angular dependence and 
a spherical Bessel function ji(m n r) providing the radial 
one {r,d,ip - spherical co-ordinates). Here, /j,f n is the cor- 
responding decay rate with n giving roughly the number 
of zeros between center and surface of the PNS. 

Table [21 summarizes the critical magnetic Reynolds 
numbers Rm crlt for the different velocities and field types 
and some characteristic properties of the marginal mag- 
netic eigensolutions, which are all non-oscillatory, that is, 
3{A} = 0. 

In Figs. |21 and we present these magnetic fields on 
the surface of the PNS and on a spherical surface in- 
side the convective shell where field generation is sup- 
posed to be most vigorous. Note that the Al solutions 
contain quadrupoles, but those without any axis of sym- 
metry (>2 (#, We emphasize that the magnetic fields 
shown belong to velocity fields extracted from the evo- 
lution of the PNS convection and then considered fixed, 
i.e., the fields corresponding to the later moment has not 
evolved from the one corresponding to the earlier. 



In their spherical harmonic spectra (referring to lati- 
tudinal scales), the 'late' magnetic fields show somewhat 
larger scales than the 'early' ones, but smaller scales in 
their spherical Bessel spectra (referring to radial scales). 
Both SI fields are at the surface predominant around the 
equator where the early field has its maxima while the 
maxima of the late one occur at ±30°. 

3. Conclusions 

The main conclusion to be drawn is that the geometrical 
structures of the velocity fields employed are well suited 
to amplify a magnetic seed field. However, given the huge 
elec tric conductivity of the PNS matter (w 10 24 s -1 , see, 
e.g.. lThompson fc Duncan! l)l993|h Table 1) and convective 
velocities of « 10 8 cm s _1 , the magnetic Reynolds numbers 
turn out to be by about 16 orders of magnitude overcriti- 
cal. As in general there is no monotonous relation between 
the growth rate of a dynamo and Rm nothing can thus 
be inferred with certainty about the dynamo activity of 
the considered convective flows at their given strengths. 
However, it is known that due to increasing flux expul- 
sion the growth rate as a function of Rm will typically 
reach a maximum and then decrease until the dynamo- 
ability of the flow may get lost completely or at least the 
growth rate appro aches zero ('slow dynamo'). [See, e.g., 
iFuchs et al ] lll997h for the Gubbins flow, where the ratio 
of the upper to the lower critical Rm for the m = 1 modes 
is only 5, or iRadler et alJ l)2002l) where it is shown that 
the a coefficient of the mean-field model of the Karlsruhe 
dynamo experiment decreases with growing fluid velocity.] 

There is no way to calculate magnetic field growth 
rates with the velocity amplitudes given above, as all ex- 
perience shows that the field scales get smaller and smaller 
with growing Rm. In our case, we are sure that with the 
available computing resources the necessary spatial res- 
olution is by far unreachable. Even if the problem would 
perhaps be less difficult in ideal MHD, one has to be aware 
that then the Bondi-Gold theorem prohibits a kinematic 
dynamo solution unequal to zero outside the PNS. 

On the other hand, PNS convection is a transient phe- 
nomenon starting from small amplitudes and ending with 
the fluid at rest. Let us discuss the possible scenarios as- 
suming that the convection starts with a sufficiently small 
Rm and then passes, continuously growing, through a 
stage which is critical with respect to the dynamo. 

First we assume that all stages of the convection be- 
yond that point and up to a second critical stage only 
close to the very end of the convective phase show posi- 
tive growth rates of the corresponding kinematic dynamos. 
Then, a magnetic field will be generated anyway (see 
Fig. |2 case (i)) 3 . 

In a contrasting scenario, the increasing convection 
could leave its dynamo-active phase very soon after hav- 

3 We disregard here the somewhat exotic possibility of a 'sui- 
cidal' dynamo due to the existence of two or more stable con- 
vective states which differ in their critical magnetic Reynolds 
numbers. Then, in principle, the magnetic field co uld switch an 
over- into an undercritical state 'killing' itself, see lFuchs et alJ 
il999Tl . 
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Fig. 2. Marginal dynamo fields for the velocity from the model of lKeill at 30 ms after bounce. Left: SI, right Al 

solution. Upper row: r — R, lower row r = 0.5R. Arrows: tangential components, color encoding: normal component. 
The field strengths are arbitrary as Eq. Q is linear and homogeneous. 



ing entered it (see Fig. 0] case (ii)). Then again two pos- 
sibilities exist: If the magnetic field grows so slowly that 
it cannot influence the flow significantly before it loses 
its dynamo-capability again, no strong fields will be gen- 
erated and the vigorous convection can even destroy the 
(weak) field it had begun to build up before (see Fig. 0] 
case (iia)). In the opposite case, the magnetic field grows 
so fast that its Lorentz forces become able to hinder the 
convection from reaching amplitudes too high for dynamo 
action (see Fig.0J case (iib)). That is, the generated field 
could fundamentally change the flow in a very early stage. 

In the first case of the second scenario, a further op- 
tion is the recovery of the dynamo ability of the convection 
during its dying-out. Then, the period during which the 
growth rate of the field is positive will be very short, limit- 
ing the accessible magnetic field strengths already on the 
kinematic level (see Fig. 0] case (iib)). 

When accepting that the second scenario — the con- 
vection is a 'slow dynamo' — is the more probable one as 
'fast dynamos' (growth rate positive and not approaching 
zero for Rm — > 00) are rarely found, we must either ques- 
tion the existence of an efficient PNS dynamo or declare 
that the influence of the generated magnetic field must 
not be neglected. Consequently, the possible essential role 
of PNS convection in enabling the supernova explosion 
must then be revised given its significant modification by 



the magnetic field at a very early stage. The only way to 
decide between the expounded alternatives is to perform 
MHD -calculations starting from the very beginning of the 
convection, i.e., with extremely small velocities. 

We concede, that the two major scenarios (cases (i) 
and (ii) in Fig. 0J| represent only the extrema of an in- 
finitude of possibilities which could be constructed when 
assuming rapid changes in the dynamo-related properties 
of the convection during its lifetime: Already for a fixed 
flow geometry the dynamo growth rate p as a function of 
Rm may change its sign several times in the relevant range 
of values. But as the flow geometry is not fixed during the 
lifetime of the convection, even its characterization as a 
slow or fast dynamo could be different at different mo- 
ments. Nevertheless, for any time interval with a positive 
p terminated by a sign change of p the discussion summa- 
rized as cases (iia) and (iib) can be applied likewise leading 
again to the above crucial alternative. Of course, possible 
consequences for the role of PNS convection in supporting 
SN explosions can only be expected during (roughly) the 
first second of the PNS's life. 
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Fig. 4. Schematic representation of possible PNS dynamo scenarios. Rm - magnetic Reynolds number (solid and 
dotted [lower right panel] /green); p = 5R{A} - dynamo growth rate (dashed/red); E mag - energy of the entire PNS 
magnetic field (dash-dot /blue). Colors refer to the online version. Note that the curves for Rm, p and E mag are sketched 
as if no back-reaction of the magnetic field on the convection existed (exception: dotted curve in lower right panel) . 



